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Abstract 

Generally, in the financial literature, the notion of quadratic VaR is implicitly confused with the 
Delta-Gamma VaR, because more authors dealt with portfolios that contained derivatives instruments. 
In this paper, we postpone to estimate both the expected shortfall and Value-at-Risk of a quadratic 
portfolio of securities (i.e equities) without the Delta and Gamma Greeks, when the joint log-returns 
changes with multivariate elliptic distribution. To illustrate our method, we give special attention to 
mixture of normal distributions, and mixture of Student t-distributions. 
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1 Introduction 

Value-at-Risk is a market risk management tool that permit to measure the maximum loss of the portfolio 
with certain confidence probability 1 — a, over a certain time horizon such as one day. Formally , if the price 
of portfolio's P(t, S(t)) at time t is a random variable where S(t) represents a vector of risk factors at time 
t, then VaR be implicitly given by the formula 

Prob{-P{t, S(t)) + P(0, 5(0)) > VaR a } = a. 

Generally, to estimate the VaR a for portfolios depending non-linearly on the return , or portfolios of 
non-normally distributed assets, one turns to Monte Carlo methods. Monte Carlo methodology has the 
obvious advantage of being almost universally applicable, but has the disadvantage of being much slower 
than comparable parametric methods, when the latter are available. 

In this paper we are concerned with the numerical estimation of the losses that the portfolio of equities 
faces due to the market, as a function of the future values of S. Following the quadratic Delta-Gamma Port- 
folio, we introduce the notion of quadratic portfolio of equities due to the analytic approximation of Taylor 
in 2 nd order of log-returns for very small variations of time. Quadratic approximations have also been the 
subject of a number of papers dedicated to numerical computations for VaR ( but these have been done 
for portfolio that contains derivatives instruments). We refer the reader to Cardenas and al.(1997) |2] for a 
numerical method to compute quadratic VaR using fast Fourier transform . Note that in P^, Brummelhuis, 
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Cordoba, Quintanilla and Seco have dealt with the similar problem, but their work have been done analyt- 
ically for Delta-Gamma Portfolio VaR when the joint underlying risk factors follow a normal distribution 
which is a particular case of elliptic distribution . All our calculus will be done according to the assumption 
that the joint securities (i.e equities) log-returns follow an elliptic distribution. To illustrate our method, we 
will take some examples of elliptic distributions as mixture of multivariate t-student and mixture of normal 
distributions. Note also that, Following RiskMetrics, Sadefo-Kamdem JT](2003) have generalized the notion 
of A-normal VaR by introducing the notion of A-Elliptic VaR , with special attention to A-Student VaR , 
but this concerned the linear portfolio. In this paper, we will do the same for nonlinear quadratic portfolios 
without derivatives instruments. 



The rest of the paper is organized as follows: In section 2, we introduced the notion of quadratic portfolio 
of securities(i.e equities) due to the 2 nd order Taylor approximation of log returns . Our calculus is made 
with the more generalized assumptions that the joint underlying log-returns follow an elliptic distribution. 
That is why in section 3, following [Ij, we recall the definition of elliptic distribution and we show that under 
the hypothesis of elliptic distribution the VaR estimation of such portfolios is reduced to a multiple integral 
equation. Next following the paper of Alan Genz f| , we recall the notions of symmetric interpolar rules 
for multiple integrals over hypersphere and we use this method to reduce our problem to one dimensional 
integral equation. In the same section, we illustrate our method by giving an explicit equation with solution 
VaR ( Value-at-Risk ) when the joint log-returns follow some particular mixture of elliptic distributions 
named mixture of multivariate Student ^-distributions or the mixture of normal distributions, in these cases 
the VaR estimation is reduced to finding the zero's of a certain specials functions. In section 5 we treats the 
expected shortfall for general elliptic quadratic portfolios of securities without derivatives instruments and 
we illustrate with the special case of normal distribution. Finally, in section 6 we give a conclusion. 



2 Quadratic Portfolio of Securities(i.e Equities) 

A portfolio of n securities is a vector 9 £ M. n ; the component 9i represents the number of holdings of the i th 
instruments, which in practise does not need to be an integer. So at time t the price of the portfolio of n 
securities is given by: 

n 

P(i) (1) 
»=1 

where S(t) = (Si(t), S n {t)) such that 

n n 

p{t) - p(o) = £ di(Si(t) - $(o)) = ^°)-^s7oT ~ 1} 

For small fluctuations of time and market, we assume that log-return is given by : 

log(S i (t)/S i (0))=r h (t) (2) 

therefore 

S z (t) - Si(0) = ^(0)(||| - 1) = Sm(expMt)) 1) 

Then we have that 

S(t) = (Si(0)exp(T}i) t S n (0)exp(ri n )). 

By using Taylor's expansion of the exponential's due to small fluctuations of returns with a time, we have 
that : 

exp(m(t)) - i « m(t) + ^4r-- (3) 
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If we assume that r\ = (771 , . . . ,rj n ) is an elliptic distribution note by N n (fi, £, <j>) then 
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n n , , 

P(t) P(0) = ^ Si(0) • 0i • (espfo,(*)) - 1) » XI Si(o)Mm(t) + v Af). (4) 

1=1 i=l 

By following the usual convention of recording portfolio losses by negative numbers, but stating the 
Value-at-Risk as a positive quantity of money, The VaR a at confidence level of 1 — a is given by solution of 
the following equation: 

Prob{\P(t) - P(0)\ > VaR a } = a 
In the probability space of losses \P(t) — P(0) \ = —P(t) + P(0), therefore recalling we have 
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Prob{J2 Si(p).0i(iH(t) + ^Y~) < -VaR a } = a, (5) 
then an elementary mathematical tool give that : 



m(t) + A = l {rii{ t) + if-i ( 6 ) 



t 

2 2 V 
therefore 

n „ n „ 

Prob{J2 SM-jMt) + I) 2 ) < -VaR a + E ^ ' S ^ = a - ^ 

i=l i=l 

By posing X — (rji +1, rj n + 1), it is straightforward that , X is an elliptic distribution due to the fact that 
it is a linear combination of elliptic distribution r\. We note X ~ N(fi + 1, S, (j) ) with a continuous density 
function /ii (x) . Remark that 

n 

^a,fe(f) + l) 2 - (x,A.a?) 

with A = (aji) i=1 „ is a diagonal matrix with diagonal values a« = Sl ^ £ * > and fi = + . . . , /x n +l) = 
fj, + I is the mean vector of X , therefore J7J) becomes 

Pro&{ pf,A.X) < k} = a (8) 

with fc = —VaR a + Y^i=i a a — ~ VaRa- We will suppose that k > 0, this means that the Value-at- 
Risk of our portfolio's is not greater than P(0)/2 . 

Remark 2.1 We remark that to estimate the Value-at-Risk of a portfolio of securities (i.e equities), the 
computation of our model need as inputs the quantity 0i and the initial security price Si(0) for each i — l..n 
as given in £7J). Recall that in literature, the computation of Quadratic Delta-Gamma (A-T) VaR need as 
inputs the sensitivity vector A and the sensitivity matrix T, because of the presence of derivatives products 
in the portfolio. 



3 Reduction to an Integral Equation 

In this section, we will reduce the problem of computation of the Value at Risk for quadratic portfolio of 
equities to the study of the asymptotic behavior of the density function distribution over the hyper-sphere. 

3.1 Notions of Elliptic Distributions 

The following definitions will be given as in 0](2OO2) . 
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3.1.1 Spherical Distribution 

Definition 3.1 A random vector X = (Xi, X 2 , X n ) 1 has a spherical distribution if for every orthogonal 
map U € i? nx ™ (i.e. maps satisfying UU 1 = WU = I n xn) 



we note that: X ~ S n (<fi).If X has a density f(x) then this is equivalent to f{x) = g{x t x) = g(\\x\\ 2 ) for some 
function g : R + — ► R + , so that the spherical distributions are best interpreted as those distributions whose 
density is constant on spheres. 

Elliptical distributions extend the multivariate normal N n (fi, E), for which /j, is mean and E is the 
covariance matrix. Mathematically, they are the affine maps of spherical distributions in M™. 

3.1.2 Elliptic Distribution 

Definition 3.2 Let T : R n — > R n , y i — > Ay + fj, , A e R nxn , \i G R n . X has an elliptical distribution 
if X = T ( Y ) and Y ~ S n (cp). If Y has a density f(y) = g(y t y) and if A is regular (det(A) ^ so that 
E = A f A is strictly positive) , then X = AY + fi has a density 



and the contours of equal density are now ellipsoids. An elliptical distribution is fully described by its mean, 
its covariance matrix and its characteristic generator. 

• Any linear combination of an elliptically distributed random vector is also elliptical with the same 
characteristic generator <j> . If Y ~ N n (p, E, , b € R m and B € R mx " then B.Y + b - N m {Bfi + 
b,BY,B\<f>) . 

3.2 Integral Equation with solution VaR 

Since X is an elliptic distribution, its density take the following form: 



UX = d X. 1 



h(x) = g((x - tfX-^x - ^))/^de~t^) 



(9) 



/iiO) = h((x-T)) 

II is the vector of unities and h is the density function of rj which take the form : 



h(x) = g((x - /iJE-^a: - p)*) / y/ det(V) , 



therefore we have the following equation 



n 



Prob{{X, AX) > -VaR a + ^ a u } = 1 



a = I(k) 



l 



with VaR a , as solution such that k = —VaR a + J2i a a- 

In terms of our elliptic distribution parameters we have to solve the following equation: 




(10) 



l{k) = I g((y m - l^s-^y - m - 

J{(x,A.x)>k} ^det(E) 



1 =d denote equality in distribution 
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since A is a diagonal matrix with all positive diagonal values, we decompose A = A 1 / 2 • A 1 / 2 therefore the 
equation (|1L)|) becomes 

i(k)= [ g(\\z\\ 2 )dz= [ g{\\z\\ 2 )dz. 

Cholesky decomposition states that £ = AA l when E is suppose to be positive, therefore if we changing the 
variable z = A~ l (y — fx — J), the precedent integral becomes : 

I(k) = [ g(\\z\\ 2 )dz. 

J (Az+^t+I) t A(Az+/i+I)>k} 

If we do the following decomposition (Az +fi+ I)*A(Az + fj, + 1) = (z + v)*D(z + v) + 6, with D = A* ■ A • A 
, v = A^ 1 (^i + I) and S = 0, after some elementary calculus 

I(k) = I g(\\z\\ 2 )dz. 

J {(z+vy D(z+v)>k-S} 

If we suggesting ki = k — S = R 2 , z + v = u, dz = du, we find that 

I{k) = I g(\\u - v\\ 2 )dz = l-a. 

J^.D.u^R 2 } 

By introducing the variable z = D 1 / 2 u/R, we have that: 

2 dz 



I{R) = R n / g(\\RD— z - v\\ 2 ) 



{|MI>i} y/det{D) 

next, by using spherical variables z = r.£ with £ S S n -i and dz = r n ~ 1 drda(t;), where da(z) is a 
elementary surface of z on S n -\ — € £i + £l + — + £n = 1} an d by introducing the function J(r, i?) 
such that 

J(r,R)= [ g(\\rRD^^-v\\ 2 )d<j(^ (11) 



we obtain 



R- n -I(R)= r ,w / ffdlri^D— C - v\\ J )da(0 , = r"~V(r,fl) (12) 



Next By introducing the function 



our goal will be to solve the following equation 



dr f°° n _ x „ dr 



y/det{D) Ji y/det{D) 



H(s) = s n / ^-^(r.ajdr, (13) 



#(s) = (1 - a)yJdet(D). (14) 

In the following section, we propose to approximate J(r,R) by applied the numerical methods giving in 
the paper of Alan Genz (2003), ( see [H] for more details). 

3.3 Numerical approximation of J(r,R) 

In this section, we estimate the integral J(r,R) by a numerical methods given by Alan Genz in [S] . 
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3.4 Some interpolation rules on S n ~i 

The paper |H| of Alan Genz, give the following method. Suppose that we need to estimate the following 
integral 

J(f) = f f(z)da(z) 

where da(z) is an element of surface on S n -i — {z\z 6 R n , z\ + z\ + ... + z% = 1}. 

In effect, let be the n-1 simplex by T n _i = {x\x 6 K™ _1 ,0 < X\ + x 2 + ■■■ + x n -\ < 1} and for any 
x G T n _i, define x n = 1 — Yn=i Xi ■ Also ~ (*pn •••>^pn-i) ^ points t , t\,. . . , f m are given, satisfying 
the condition : \t p \ = X)"=i ^p; = 1 whenever "Y^l—xPi = m i f° r non-negative integers pi,. . . , p n , then the 
Lagrange interpolation formula (sylvester |12j for a function g(x) on T n _i is given by 

n Pi-1 2 _ f 2 

l^( 9>X )= Enn^ 

|p|=mt=l i=0 Pl -J 

L( m ' n_1 )(p,x) is the unique polynomial of degree m which interpolates g(x) at all of the C^ +n _ 1 points in 
the set {x\x — (t pi , t Pnl ), \p\ = to}. Silvester provided families of points, satisfying the condition \t p \ = 1 
when \p\ = to, in the form t{ = ^+8^ for i=0,l,. . . ,m, and /i real. If < 6 < 1 , all interpolation points 
for L( m > n - 1 )(s,x) are in T n —\. Sylvester derived families of interpolators rules for integration over T„_i 
by integrating L^ m,n ^ 1 \g,x) . Fully symmetric interpolar integration rules can be obtained by substitute 
Xi = zf , and ti — uf in L^- m ' n ~ 1 \g, x), and define 

npi — l 2 2 
|p|=m i=l j=0 ^ 

where f{u} is a symmetric sum defined by 

f{u} = 2~ c(tl) y^J{sxUx,s 2 U2, ■ ■ .,s n u n ) 

s 

with c(u) the number of nonzero entries in (u\, . . . , u n ), and the ^ s taken over all of the signs combinations 
that occur when Si = ±1 for those i with different to zero. 

Lemma 3.3 If 

n Pi — 1 2 _ 2 

i=l j=Q P^ J 

then 

J{f)=R (m,n) {f)= J- Wp f{u p } 

\p\—m 

f{u} = 2~ c(u) ^2 f(siui,s 2 u 2 , s n u n ) 

s 

, with c(u) the number of nonzero entries in (u\, ...,«„), and the ^ s taken over all of the signs combinations 
that occur when Sj = ±1 for those i with ui different to zero. 

The proof is given in [H] by (Alan Genz (2003)) as follow: 
Let z k — z kl Zj 2 • z* n . J and R are both linear functionals, so it is sufficient to show that t& n < n \z k ) — J(z k ) 
whenever |fc| < 2to +1. If k has any component ki that is odd, then J(z k ) = 0, and R^ m ' n ^(z k ) — because 
ever term u k in each of the symmetry sums f{u p } has a cancelling term - u k . Therefore, the only monomi- 
als that need to be considered are of the form z 2k , with \k\ < to. The uniqueness of L^ m,n ~ 1 \g,x) implies 
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l\m,n ; x ") _ x k w henever |fc| < m, so M^" 1 -™) (z 2k , z) — z 2k , whenever |fe| < m. Combining these results: 

J(/)=M^ (/,«))= £ w p f{u p } = R^ n \z k ) 

\p\=m 

whenever f(z) = z k , with |fc| < 2m + 1, so R^ m ' n \f) has polynomial degree 2m+l. For more details ( cf. 
Genz 0). 

3.5 Application to Numerical approximation of J(r, R) 

Since our goal is to estimate the integral (|10fl . it is straightforward that the theorem (4.2.1) is applicable to 
the function f such that 

f(z)=g(\\rRD^z-v\\ 2 ). 

Next, since we have that 

f{u p } = g(\\rRD^( S .u P Y ~v\\ 2 ), 
by introducing the approximate function J Up that depend to the choice of u p , H1(J|) becomes 

J(r,R)^ Yl J2 w p giWrRD^is.UpY - v\\ 2 ) = J Up (r, R) (15) 

we note s.u p = (sxUx, s n ii n ) . 

Remark 3.4 J u (r,R) is the numerical approximation of J(r,R) as given in f!5fjl. is depend to the choice 
of interpolation points u p on hypersphere. Recall that J(r, R) was a fixed function that depend to R and the 
density function of our elliptic distribution . 

By introducing H Up , the approximate function of H as define in that depend of J Up , such that 

/oc 
r n - 1 J Up (r,s)dr^H(s). (16) 

By replace H(s) in (|14|) by H Up (s) we then prove the following result: 

Theorem 3.5 If we have a quadratic portfolio of securities (i.e equities) such that the Profit & Loss function 

over the time window of interest is, to good approximation, given by All « ' ^i{Vi{t) + )' 

with portfolio weights 0i. Suppose moreover that the joint log-returns is a random vector (771, . . . ,r) n ) that 
follows a continuous elliptic distribution, with probability density as in @, where fi is the vector mean and 
S is the variance- covariance matrix, and where we suppose that g(s 2 ) is integrable over M, continuous and 
nowhere 0. Then the approximate portfolio 's quadratic elliptic VaR^ Up at confidence (I — a) is given by 

VaRZ tUp = ^p--Rl Up (17) 



where Rg lU „ is the unique solution of the equation 



(I -a) 



H Up (s) = (!-«)• y/dei(Dj = ^—^ 



\ 



det(T,)J[ei ■ Sj(0). (18) 



i=i 



In this case, we assume that our losses will not be greater than half-price of the portfolio at time 0. 

Remark 3.6 The precedent theorem give to us an approximate Quadratic Portfolio Value-at-Risk (V aR^ u ) 
that depend to our choice of interpolation points on hypersphere ,a and the function g. Therefore it is clear 
that the best choice of interpolation point will depend to the g function in f^) . 
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With some simple calculus we have the following remark 
Remark 3.7 

J Up (r,R)= ^2 w p ^g{a(s,u p ,R)-r 2 -2-b{s,u p ,R,D,v)-r + c{v)^j (19) 



\p\ — m 



for which a(s,u p , R, D) = \\RD ? (s.u p )\\ 2 , b(s,u p , R, D,v) = R < D 2 (s.u p ), v >,c = \\v\\ 2 . Sometimes, for 
more simplification we will note a,b,c. 

Since inequality of Schwartz give that b 2 — ac < , we use the change of variable by posing 61 = b ~ ac < 0, 
u = r — by using using the binom of Newton, and by introducing the function G*' 9 for j = 0, ..,n — 1, 



such that we have the following remark 
Remark 3.8 



71 — 1 / - 

n — 1 



J(r, R)=Y, ™pEE ■ (b/af^G^JR) (20) 

|p|=m s j=0 V / 
/■oo 

G^JR)= z> ■g(az 2 -b 1 )dz (21) 



/or which a,b and c are defined in \3. 7| ). 

-i 

By replace d = b/a = — 2 _p-"p)- t ' > by its value in 13.81) .we obtain 

i?i|DT-(s.u p )|| 2 



G%,{R)= <D ^ ( . > ^^ g (i? 2 ^|^(^ P )||V- < ; n i S ; M ^ 2 > + lM| a )<fe (22) 
■/l- — 2 _ l ' a "p'- t,> V ||£)— (s.itp)|| 2 7 

R||DT ( s . Up )|| 2 

we then have the following theorem 



Theorem 3.9 If we have a quadratic portfolio of securities (i.e equities) for which the joint securities log- 
returns changes with continuous elliptic distribution with pdf distribution as in ifffil. then the approximate 
portfolio 's quadratic elliptic VaR 9 a Up at confidence (1 — a) is given by 

VaRi <Up = - Rl Up (23) 
where R g , u is the unique solution of the equation 



n— 1 - 1 n—l—j 

E E E n 7 • (S^?) • - — 



det(E)Jp< •£<(()). (24) 



In this case, we assume that our losses will not be greater than half-price of the portfolio at time 0. 

Remark 3.10 We have reduced our problem to one dimensional integral equation. Therefore, to get an 
explicit equation to solve, we need to estimate G^ 9 S (R) that depend to R with parameters g,u p ,v and D. 

Therefore, in the case of normal distribution or i-distribution, it will suffices to replace g in the expression 
of H21(l . an to estimate the one dimensional integral J2U- 
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3.5.1 The case of normal distribution 

In the case of normal distribution, the pdf is given by: 



fix) = —=J- exp(-\(x - ulirVx - fif) (25) 



and specific is given as follow 

g(x) = {2ir)~^e~^ = C(n).e" 
therefore it suffices to replace g in <|2"T|) then 



Gi p . s (i?) = (27r)-^e^ / u=e~^du. (26) 
If 1 — — > ( it is the case when R is sufficiently big such that \v\ < R\\D~ (s.u p )\\ ). 



G kA R ) = ( < D^js.up), v > 2 i+i j + 1 (fl||J^(aS)||) 2 <D^(s.u p ),v> 

ex p(-WP)(27r)-f E||D^(s%)|| 2 2 ' 2 fl ||£,^ (s^H 2 

(27) 

therefore, since a = (R\\D^~ (s.u p )\\) 2 we have the following theorem 

Theorem 3.11 If we have a portfolio of securities (i.e equities), such that the Profit & Loss function over 

the time window of interest is, to good approximation, given by All w 2i=i Si(0) ■ Oi(r]i(t) + 2 ' ), with 
portfolio weights Oi. Suppose moreover that the joint log-returns is a random vector (771, ... , r) n ) that follows 
a continuous multivariate normal distribution with density function in 125\) . vector mean [i , and variance- 
covariance matrix X, the Quadratic Value-at-Risk (VaR a ^ Up at confidence 1 — a is given by the following 
formula 

2 P (°) 
R Up , a - ~VaRa,u p H — 

for which R Up ,a is the unique solution of the following transcendental equation. 

for which h = «^^/ s -"p)."» 2 _ |i w ||2 b = <D-rjs.u p ),v> = R 2\\ D ^r s . u )f. /„ ftfc case, we implicitly 
assume VaR a:Up < P(0)/2. T is £/ie incomplete gamma function. 

3.5.2 Case of i-student distribution 

If our elliptic distribution is in particular chosen as the multivariate t-student distribution, we will have 
density function given by 

«w = 1 + -) w - c(„,„)(i + £) w (2 9) 



r(i//2).7T"/ 2 V 

therefore by replacing g in i|24[l . we obtain the equation 



n-l 



|p|=m s j=0 
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suggesting c\ = v — b\ , our equation is reduce to 

E ^EEMw^ 4 ^(au^ + c 1 ) i ^ ) du= i -^hlM£l. (31) 



\p\=m s j=0 



3 J Ji-i v ' v—C{n,u)R n 



changing variable in this integral according to v = u 2 and (3 = , we find that 



(32) 

For the latter integral equation, we will use the following formula from [7]: 
Lemma 3.12 (cf. 'V, formula 3.194(2)). If \arg{f)\ < n, and Re(v{) > Re(fi) > , then 

+oo u^ Vl B~ Vl 1 

x^- 1 (l + fSxy^dx = '- 2^1(^1,^1 - fJ.;vi -n+1;-- — ). (33) 

V\ — fl p ■ u 

Here 2^1 (ck; f3, 7; w) is the hypergeometric function. 

In our case, 2^1 = !i P^, u = (1 — b/a) 2 , V\ — fi = SdJZzink^ — [i + \ = n+lJ ~'J +1 therefore If we replace in 
(|32[1 . we will obtain the following result. 

Theorem 3.13 If we have a portfolio of securities (i.e equities), such that the Profit & Loss function over 

the time window of interest is, to good approximation, given by ATI « S"=i ' @i(Vi(t) + ^ ^ ) , luii/i 

portfolio weights 6i. Suppose moreover that the joint log-returns is a random vector (771, ... , r) n ) that follows 
a continuous multivariate t- distribution with density function given by vector mean \jl , and variance- 

covariance matrix £, the Quadratic Value-at-Risk {VaR a , Up at confidence 1 — a is given by the following 
formula 

R 2 VnR + P(0) 

for which R Up ,a is the unique solution of the following transcendental equation. 

1 77' n+v n+v—j — 1 . n+v—j + l . 6? — v 

R n v yyM (V^ 2 M~ — s — ■ — 2 — _ (7r/2)"/ 2 rW2) 



ft- a )|pNm , ^ W (1/ - 6?)^ (n + * ~ 3 ~ 1 Vdet(E) EHU ^ • 5*(0) ^r(^) 

— 1 2 — 1 

for which bi = ^ <D 2 _ 1 ( s -"p)' t,> ) ||u|| 2 , - = — 2 _p-"p)' , " > a = R 2 \\D~ (s.u v )\\ 2 . In this case, we implicitly 

\\D^r{s.u P W a R\\D^r(s.u P W 

assume that VaR aAlp < P(0)/2 

Remark 3.14 Note that, Hypergeometric 2-Fi's have been extensively studies, and numerical software for 
their evaluation is available in Maple and in Mathematica. 



4 Quadratic VaR with mixture of elliptic Distributions 

Mixture distributions can be used to model situations where the data can be viewed as arising from two 
or more distinct classes of populations; see also [§]. For example, in the context of Risk Management, if 
we divide trading days into two sets, quiet days and hectic days, a mixture model will be based on the 
fact that returns are moderate on quiet days, but can be unusually large or small on hectic days. Practical 
applications of mixture models to compute VaR can be found in Zangari (1996), who uses a mixture normal 
to incorporate fat tails in VaR estimation. In this section, we sketch how to generalize the preceding section 
to the situation where the joint log-returns follow a mixture of elliptic distributions, that is, a convex linear 
combination of elliptic distributions. 
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Definition 4.1 We say that {X\, . . . ,X n ) has a joint distribution that is the mixture of q elliptic distri- 
butions N(^ij,T,j,(j>j) 2 , with weights {f3j} (j=l,..,q ; f3j > ; 5Z? =1 j8j = -0' ^ ^ s cumulative distribution 
function can be written as 

<? 

Fxi,...,X n (%i, ...,x n ) = ^2f3 Fj{xi, . . . ,x n ) 
j=i 

with Fj{x\, . . . , x n ) the cdf of N(jij, Ej, (^j). 

Remark 4.2 In practice, one would usually limit oneself to q = 2, due to estimation and identification 
problems; see [5]. 

We will suppose that all our elliptic distributions N(p,j, Sj, <j>j) admit a pdf : 

f 3 {x) = ^-^g^x-^f^x-^f) (35) 

for which each gj is continuous integrable function over R, and that the gj never vanish jointly in a point of 
K 9 . The pdf of the mixture will then simply be Ylj=i Pjfj( x )- 

Let 

= Aj Aj 

So, following l(l"2")l.we introduce J k (r,R) such that 

air " = E/ ^ 9k (\\rRD^Z-v k \\ 2 )d*(0 ^===Y, r n ~\j k {r,R)dr (36) 

fc=] yi L "'S„_ 1 J \Jdet{D k ) k=1 Ji 

Next following Q21[l. we introduce the function 

GZs,ki R )= Rn b z ' -gk(a k z 2 -b lk )dz (37) 
j i — 

_! b 2 —a c 

witha fe = \\RD—(s.u pk )\\ 2 , b k = R < D k 2 (s.u pk ),v >,c k = \\v k \\ 2 ,bi k = k a ° fcCfc , then we have the following 
corollary 

Theorem 4.3 If we have a portfolio of securities (i.e equities) such that the Profit & Loss function over 

the time window of interest is, to good approximation, given by All X)"=i <%(0) ' @i(Vi(t) + 2 ) ' w ^ 
portfolio weights Qi. Suppose moreover that the joint log-returns is a random vector (771 , . . . ,T) n ) is a mixture 
of q elliptic distributions, with density 

Kx) =J2(3 j \X j \- 1/2 g j ((x - n)Yq\x - h)*) 

where [ij is the vector mean, and Ej the variance- covariance matrix of the j-th component of the mixture. 
We suppose that each gj is integrable function over K, and that the gj never vanish jointly in a point of 1R™ 1 . 
Then the value- at- Risk, or Quadratic mixture-elliptic VaR ; at confidence 1 — a is given as the solution of 
the transcendental equation 



71-1 



k=l \p\=m s j=0 



n-l\ rh/ , m _ 1 - J -^. fc (ff-^^) 1/2 )_(l 



y y (w"' hj!! p „ — _^ = i±_^ (38) 



/or which G J ^ 3 s k is defined in \31\j . In this case, we assume that our losses will not be greater than half-price 
of the portfolio at time 0. 

2 or N(fij , E j , gj) if we parameterize elliptical distributions using gj instead of <pj 
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Remark 4.4 One might, in certain situations, try to model with a mixture of elliptic distributions which 
all have the same variance-covariance and the same mean, and obtain for example a mixture of different tail 
behaviors by playing with the g^s. 

The preceding can immediately be specialized to a mixture of normal distributions: the details will be 
left to the reader. 



4.1 Application with mixture of Student i-Distributions 

We will consider a mixture of q Student ^distributions such that, the k th density function i — 1, .., q will be 
given by 

»<»> - jm„(l + - C M {l + •)« = *«> (39, 

T(v k /2).TT n l Z \ V k ) \ V k J 

and Sfc = A\A k therefore by replacing g k by g in H4()(l and since integration is a linear operation, we obtain 
the following theorem 

Theorem 4.5 If we have a portfolio of securities (i.e equities) such that the Profit & Loss function over 

the time window of interest is, to good approximation, given by All rs Y^j—i Si(0) ' + ^2 )> w ^ 

portfolio weights 9i. Suppose moreover that the joint log-returns is a random vector (rji, . . . , 7?„) is a mixture 
of q t- distributions, with density 

« . T(^) / (x- uA^J^x- aj)\-^p- 

- g^i^i 7 r^Tsf^C 1 + ; — ) 

where \ij is the vector mean, and Ej the variance-covariance matrix of the j-th component of the mixture. 
We suppose that each gj is integrable function over K, and that the gj never vanish jointly in a point of W 71 . 
Then the value- at-Risk, or Quadratic mixture-student VaR, at confidence 1 — a is given by : 

< = ~VaR a + ^ 
for which R Up is the unique positive solution of the following equation: 



WpT (^±R) n ^ f n -l\ R n( bk / ak )n-i-j 2Fl 



n+v k n+Uk—j—1 . n+u k -j + l . b lk -b 

2 ' 2 ' 2 ' rrm 



Q fc( 1 --^f 



h_\2 



(1 



h\tm r( ^ /2) - sh^ j J ((^-bi k yu k )^ (n+^-j-i)vmn^AV) - (7r/2 ) 

(40) 

for which b lk = « D ~^)^» 2 _ || Ufe ||2 ) £ = <D-^s.u pkhVk > ^ ^ = ^^^^ md def ^ k) = 

\\D k 2 (s.u pk )\\ 2 R\\D k 2 (s.u ph )\\ 2 

|Sfe| . In this case, we implicitly assume that our losses will not be greater than P(0)/2. 



5 Elliptic Quadratic Expected Shortfall for portfolio of securities 

Expected shortfall is a sub-additive risk statistic that describes how large losses are on average when they 
exceed the VaR level. Expected shortfall will therefore give an indication of the size of extreme losses when 
the VaR threshold is breached. We will evaluate the expected shortfall for a quadratic portfolio of securities 
under the hypothesis of elliptically distributed risk factors. Mathematically, the expected shortfall associated 
with a given VaR is defined as: 

Expected Shortfall = E(-AII| - An > VaR), 
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see for example [§]. Assuming again a multivariate elliptic pdf f(x) — |E| 1 g{{x — /i)E x (a; — /it)*), the 
Expected Shortfall at confidence level 1 — a is given by 

-ES a = E(AIT | An < -VaR a ) 

= -E (AH • l {A n<-VaR a }) 
a 



1 

a 

|E| 



o 



((x,A.a;) -P(0)/2) ^(x) 

{(x,A.x)-P(0)/2<-VaR a } 
1/2 

((at, A.x) - P(0)/2) flr((at - /x - l)E _1 (a: - j"- l)*)da;. 

{O,A.z)<-Va.R Q +P(0)/2} 

/■Ml, 



Using the definition of VaR a and by replace All = (A, A. A) with random vector X define in section 

2, 

ES a = - ^ f (x, A.x) g((x - /x - \)TT x {x - fi - l)')dx (41) 

2 a J {{x,A.x)<-VaR, a +P(0)/2} 

Let E = A* A, as before. Doing the same linear changes of variables as in section 2 and section 3, we 
arrive at: 



ES a = 



P(0) 


R n+2 


D\ 


-1/2 


2 




a 




P(0) 


R n+2 


D\ 


-1/2 


2 




a 




P(0) 


R n+2 


D\ 


-1/2 


2 




a 




P(0) 


R n+2 


D\ 


-1/2 


2 




a 





,1 



g(\\rRD^^-v\\ A )da(0 

S„_i 



dr 



|p|=m s 

By introducing the function Q 9 Uv s such that 



r n+1 J(r, R)dr 
r n+1 J Up (r,R)dr 
V V Wp j r n + l g{\\rRD^{s.u p f-vf) dr 



(P)=P n+2 / r n+1 S (||ri2D^ -«|| 2 ) dr (42) 



we have the following theorem 

Theorem 5.1 Suppose that the portfolio is quadratic in the risk-factors X = (Ai,- • • , X n ): AIT = (A, A 
X)-?@L and that X ~ N(fi + 1, E, (f), with pdf f{x) = {E^g^x - fi - l)E _1 (a; — fJ,— If). If the VaR 
is given, then the expected Shortfall at level a is given by : 



ES a 



P(0) \D\- 1/2 _ ^ / P(0) 



2 a 

\p\—rn s 



EE™* qM (£ F-^ )1/2 )- (43) 



we introduce If and If such that 

pi pOO poo 

R- n - 2 Q 9 u P A R ) = / r n+1 g(ar 2 - 2br + c)dr = - = I^R) - /f, Up , s (P) 
Jo Jo Ji 



Following the Integral 121|) 



13 



n+1 



3=0 



for which a,b and c are defined in remark l|3.7H and 



iL p JR) = f 

Jo 



r n+ 'g{ar z - 2br + c)dr 



5.1 Expected Shortfall with normal distribution 

In the case of normal distribution, the pdf is given by l|25|l and the specific g is given as follow 

g(x) = (27r)~T e -§ = C(n).e~* 
therefore it suffices to replace g in <|2*T|) then 



Gi p , s (R) = (2 7r )-fexp(^^)y i 6 ^exp(-^)d 



■b 2 — ac\,„, .l±i + l a, 6 S , 

• S-.oC 1 -- 

2 2 a 



(2 7 r)-t CX p(^-^)(2/a) i * i r(^ 



By using the following lemma 

Lemma 5.2 (cf. formula 3462(1)). If Re(u) > 0, and Re(/3) > , then 



/ x 1 " 1 exp (-/3a; 2 - Xx)dx = {2py v/2 T{v) exp ( — ( 
Jo V8p/ V 



Here ©_„ is £/ie parabolic cylinder function with 



.„(*) = 2^V 



r(±±-) * V2' 2' 2 ; r(|) * V 2 ' 2' 2 
where $ is £/ie confluent hypergeometric function (for more details see ^ page 1018). 
we next obtain 

W» ^ S (" ; ') • «p (^) (2/«)*r(i±i, |„ - if) 



3=0 



and 



^ P , S W = (2 7 r)-texp(--^)^V«+ 1 exp(-^- T ---) 



dr 



(2tt)-9 exp(--^-)a^r(n + 2) exp (^-)b_„_ 2 (-E 



for which D_„_2 is the parabolic cylinder function. We have therefore prove the following result: 
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Theorem 5.3 Suppose that the portfolio is quadratic in the risk-factors X = (Xi, ■ • ■ , X n ): 

All = (X, A • X) - ^21 an d that X is a multivariate normal distribution. If the VaR a is given, then the 

expected Shortfall at level a is given by : 

Es a = ^- Rn+2 M~ 1/2 £ T, W P [ I UA R )- I L P ,( R )}- ( 45 ) 



\p\—m s 



for which R — \J -^j^ — VaR a 



The preceding can immediately be specialized to a mixture of normal distributions. The details will be 
left to the reader. 



5.2 Student ^-distribution Quadratic Expected Shortfall 

Following the precedent section 3,4,5, and particularly the lemma (I33|l , the application can be specialized 
to a Student t-distribution. The details will be left to the reader. 



5.3 How to choose an interpolation points u p on hypersphere 

In order to obtain a good approximation of our integral, one will choose the points of interpolation u p of our 
g function such that our approximation is the best as possible. In the case where the g function decreases 
quickly with all its derivatives of all order, in inspiring of the classical analysis, one will choose the points 
which render the maximum function \\rRD~ 1 ^ 2 (s.u p ) — v\\. 



6 Conclusion 

By following the notion of Delta-Gamma Portfolio that contains derivatives instruments, we have introduced 
a Quadratic Portfolios of securities (i.e equities) without the use of Delta and Gamma. By using the 
assumption that the joint securities log-returns follow a mixture of elliptic distributions, we have reduced 
the estimation of VaR of such quadratic portfolio, to the resolution of a multiple integral equation, that 
contain a multiple integral over hypersphere. To approximate a multiple integral over hypersphere, we 
propose to use a numerical approximation method given by Alan Genz in|S]. Therefore, the estimation of 
VaR is reduced to the resolution of one dimensional integral equation. To illustrate our method, we give 
special attention to mixture of normal distribution and mixture of multivariate t-student distribution. In 
the case of t-distribution, we need the hypergeometric special function. For given VaR, we also show how 
to estimate the expected shortfall of the Quadratic portfolio without derivatives instruments, when the risk 
factors follow an elliptic distributions and we illustrate our proposition with normal distribution by using 
the parabolic cylinder function. Note that this method will be applicable to capital allocation, if we could 
consider an institution as a portfolio of multi-lines businesses. In the sequel paper, we will dealt with this 
numerical Quadratic method, when the Portfolio contains derivatives instruments (i.e options). A concrete 
application need the estimation of w p such that \p\ — to, therefore we send the reader to Alan Genz [Sj. 
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